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ANALYTIC  TWO-DIMENSIONAL  SUBCENTER  MANIFOLDS 
FOR  SYSTEMS  WITH  AN  INTEGRAL 

Al  Kelley 

For  a  real  analytic  system  of  ordinary  differential  equa¬ 
tions  with  an  integral  H 

H  =  i  X(x*  +  y*)  +  F(z)  +  G(x,  y,  z) 
x  =  ky  +  X(x,  y,  z) 
y  =  —Xx  +  Y(x,  y,  z) 
z  =  Bz  +  Z(x,  y,  z) 

where  x  and  y  are  scalars;  z  is  an  m- vector;  X,  Y,  Z  are 
power  series  with  no  constant  or  linear  terms  \  B  is  a  constant 
matrix  with  eigenvalues  (tt,  and  iirlPi  =£  integer  (j  — 

1,  ••  •,  m);  the  existence  of  a  unique,  local,  real  analytic,  two- 
dimensional,  invariant  subeenter  manifold 

MK  =  {(*,  y,  z)  1 1  *  |  +  |  y  |  <  S,  z  =  w{x,  y )} 

is  proved,  where  w  is  a  real  analytic  function  with  no  con¬ 
stant  or  linear  terms  in  its  expansion  about  the  origin.  The 
manifold  Mx  is  composed  of  a  nested,  one-parameter  {p  S  0) 
family  of  periodic  orbits,  and  as  p-*0  the  corresponding 
periodic  orbit  goes  to  the  origin  and  its  period  goes  to  2k  |  X~l  (. 

In  this  paper  we  extend  the  result  of  C.  L.  Siegel  [7]  to  analytic 
systems  of  ordinary  differential  equations  with  an  integral  (see  (1) 
below).  C.  L.  Siegel  proved  Theorem  1  below  with  the  additional 
restrictions  that  system  (1)  is  Hamiltonian  and  that  the  matrix  B  in 
(1)  is  a  diagonal  matrix.  The  removal  of  these  two  restrictions  is  a 
technical  improvement;  our  emphasis  is  on  the  invariant  manifold 
aspect  of  the  theory  (see  Theorem  2  below).  As  far  as  proving  the 
existence  of  a  local,  nested,  one-parameter  family  of  periodic  orbits 
is  concerned,  the  theorem  proved  by  C.  L.  Siegel  [7]  is  equivalent  to 
a  result  of  Liapounov  [5]  proved  before  the  turn  of  the  century. 
However,  in  addition  C.  L.  Siegel  proved  that  the  invariant,  two- 
dimensional  manifold  comprising  the  nested,  one-parameter  family  of 
periodic  orbits  is  an  analytic  manifold,  and  that  by  means  of  an 
analytic  change  of  variables,  the  differential  equations  on  the  invariant 
manifold  have  a  very  simple  form.  Using  the  proof  method  of 
Liapounov,  one  is  able  to  conclude  only  that  the  invariant  manifold 
of  periodic  orbits  is  continuous.  (See  A.  Kelley  [3]  for  a  C*  version 
of  these  ideas.)  For  real  analytic  systems  with  an  integral  we  prove 
in  §  4  below  that  the  invariant  manifold  is  directly  computable  from 
the  differential  equations,  and  hence  unique  and  real  analytic.  Since 
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our  manifold  is  imbedded  in  the  center  manifold  (see  A.  Kelley  [4]), 
we  call  it  a  subcenter  manifold. 

The  convergence  proof  in  Theorem  1  below  is  based  on  techniques 
developed  in  C.  L.  Siegel  [7],  J.  Moser  [6],  and  A.  Kelley  [1],  [2]. 

2.  Notation.  Let  *  =  (*«  *  •  •,  z,)  be  a  complex  vector  and  let 
q  be  a  nonnegative  integer.  We  define  F«(z)  to  be  the  class  of 
functions  {F}  which  satisfy  (i),  (ii),  (iii)  below. 

(i)  F  =  F(z)  is  a  complex  scalar  valued  function,  vector  valued 
function,  or  matrix  valued  function  (dim  F  unspecified)  defined  for  z 
in  some  neighborhood  of  z  =  0. 

(ii)  F  has  a  convergent  power  series  expansion  in  z  (components 
of  z)  about  the  point  z  =  0;  the  coefficients  in  the  expansion  are 
either  scalars,  vectors,  or  matrices  (as  the  case  may  be);  there  are 
no  terms  of  degree  less  than  q  in  the  expansion  of  F. 

(iii)  For  some  K  >  0  the  components  of  F  are  majorized  by 

K(zl  +  •  •  •  +  z,)'[l  —  K(Zi  +  •  •  •  4-  z,)]-1 . 

The  constant  K  depends  on  F. 

Primes  in  this  paper  either  designate  majorants  or  designate 
variables  and  functions  in  a  complex  setting  as  contrasted  to  a  real 
setting.  Given  F(z)eF*(z)  there  exists  F’(z)  such  that 

F(z)  «  F'(z) 

where  the  symbol  <  means  F  and  F'  belong  to  the  same  r  class, 
dim  F  =  dim  F\  and  F  is  majorized  componentwise  by  F'  in  the 
sense  of  Cauchy. 

The  norm  |  ♦  |  is  the  euclidean  norm  on  vectors  and  the  operator 
norm  on  matrices.  If  F  =  F(z)  is  a  vector  valued  function  of  the 
vector  z,  then  F,  represents  the  usual  Jacobian  matrix  of  partial 
derivatives. 

If  Fe  F°(x,  y)  where  x  and  y  are  scalar  variables,  then  [FJ«(, 
represents  the  coefficient  of  the  x“yf  term  in  the  expansion  of  F. 

Let  Fe  r*(x,  y)  be  a  vector  valued,  real  analytic  function  of  the 
scalar  variables  x  and  y,  and  let 

F(x,  y)  =  £; :,mtF..,xmy'  . 

In  §4  below  we  will  use  the  notation 

Fu  =  IF ]„  =  (F.,o,  F,_,,„  -'(FJ  . 


Thus  F.  is  a  vector  in 

Throughout  this  paper  R  =  R(p,  q)  is  used  to  designate  an 
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anylytic  function,  although  occasionally  (as  above)  the  notation 
zeiJ*  is  used  to  indicate  that  z  is  an  n- vector.  The  distinction  will 
always  be  clear  from  the  context. 

3.  Analytic  changes  of  variables.  Consider  the  complex  analytic 
system  of  ordinary  differential  equations  with  an  integral  H 

H(x,  y,  z)  =  iXxy  +  F(z)  +  G(x,  y,  z ) 
i  =  iXx  +  X(x,ytz) 

V  -  -i\y  +  Y(x,y,z) 
z  =  Bz  +  Z{x,  y,  z) 

where  x  and  y  are  complex  variables;  z  is  a  complex  m- vector 
variable;  \  is  a  nonzero  complex  number;  (in  the  most  interesting 
case  X  is  real,  see  (19)  below  for  the  real  counterpart  to  system  (1)) ; 
B  is  a  constant  matrix  in  Jordan  canonical  form ;  B  =  B°  +  Bl, 

B°  =  diag(/i„  ••*,£.) 

Bl  =  subdiag  (S„  •  •  • ,  $„_,) 

with  Sj  =  0  or  d,  8  >  0  sufficiently  small  (to  be  specified  below) ; 
X,  Y,  ZeHi,  y,  z);  dim  X  =  dim  x,  etc.;  F  is  a  quadratic  form  in 
the  components  of  z  (F  &  0  is  allowed) ;  G  e  r*(x,  y,  z)  (which  implies 
that  Hert(x,y,z));  H  is  an  integral  for  system  (1). 

Theorem  1.  If  y  =£  integer  (j  —  1,  •  •  • ,  m),  then  for  system 

(1)  there  exists  a  change  of  variables  ( not  unique) 

x  =  P  +  P(P,  q) 

(2)  V  =  q  +  Q(p,  q) 

z  =  r  +  R(p,  q) 

such  that 

p  =  tX(l  +  a(pq))p  +  P(p,  q,  r) 

$  =  -»M1  +  a(PQ))Q  +  Q(P,  Q,  r) 
r  =  Br  +  R(p,  q,  r) 

P,  0,  R  =  0  when  r  =  0  , 

where  a  =  a(pq)  is  a  function  of  the  product  pq ;  a(a>)  e  /"(a))  ( which 
implies  a  {pq)  e  r*(p,  q)) ;  P,Q,Re  r*(p,  q) ;  P,Q,Re  r‘(p,  q ,  r). 

Proof.  First  we  construct  a  formal  change  of  variables  (2)  such 

that 
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i>  =  tX(l  +  a(pq))  +  P(p,  q,  r) 

-4)  4  =  +  Z(pq))  +  Q(p,  q,  r) 

1  t  =*  Br  +  R(p,q,r) 

P,  0,  &  =  0  when  r  =  0  , 

where  a  and  3  are  formal  power  series  in  the  product  pq  with  no 
constant  terms.  Then  we  will  use  the  fact  that  H  is  an  integral 
for  system  (1)  to  prove  that  3  =  a  as  formal  power  series.  Finally, 
convergence  will  be  proved  using  the  Cauchy  method  of  majorants. 

If  such  a  change  of  variables  (2)  exists,  then  from  (1),  (2),  (4) 
we  obtain 

P  =  -ixap  +  iXP  +  X(p  +  P,q  +  Q,r  +  R)  -  P{p,  q,  r) 

(5)  Ci  =  iX3q  -  iXQ  4-  Y(p  +  P,q  4-  Q,r  +  R)  -  Q(p,  q,  r) 
ft  =  BR  4-  Z(p  +  P,q  +  Q,r  +  R)  —  P(p,  q,  r) 

where  P  =  P(p,  q),  etc.  Taking  the  derivative  on  the  left  side  of  (5) 
implicity  and  using  the  condition  P,  Q,  K  ==  0  when  r  —  0,  we  obtain 


pP,  —  qP,  —  P 

=  —ap  —  apP,  +  aqP,  -  iX~lX(p  4-  P,q  +  Q,  R) 

pQ,  -qQ'  +  Q 

(6)  -aq  -  apQ,  +  3qQt  -  i\"1  Y(p  +  P,q  +  Q,  R) 

P~-j~  ~ 

=  —iX  1di_lRj^l  —  op— — Rj  4*  uq-z — Rj  iX~lZj(p  4-  P,  q  +  Q,  R) 

op  oq 

where  80  =  0.  From  (6)  the  coefficients  of  P,  Q,  R,  a,  3  can  be  com¬ 
puted  recursively.  Let 


P(P,  Q)  =  S  P*,fP'qfi  , 

a+Jkt 


etc.  From  (6)  we  have 


-  0  -  l)P«.rf>V  =  -5>a#P*V  +  ••• 

(7)  E(«  - 0  +  1)  Q..,pV  =  za..«p*fl‘+l  +  — 

£(«  —  0  +  ix~lft i)Rj,a,f  =  •  •  •  (i  =  1,  •••,»») . 


If  we  equate  corresponding  coefficients  on  the  left  and  right  side  of 
(6)  (equivalently  (7)),  we  obtain 
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(a  —  &  — 

=  l-ap  -  apP,  +  3gP,  -  i\~lX(p  +  P,q  +  Q,  £)]a,, 

(a  -  $  +  1)Q„., 

=  13 q  -  apQ,  +  SqQ,  -  i\~lY(p  +  P,q  +  Q,  P)]a>, 

(a  -  £  +  iX-'/iJR,,,', 

—  £  ap-^-Rj  +  aq-^-Rj  —  i\  1Zi(p  +  P,q+Qt  R) J  ^ . 

Clearly,  the  expressions  on  the  right  side  of  (8)  involve  only  the 
coefficients  of  P,  Q,  R  of  order  less  than  a  +  j3  and  the  coefficient  of 
Ru  •  •  • ,  Rj-i  of  order  less  than  or  equal  to  a  +  /3.  When  a  =  &  +  1, 
then  at',  must  be  chosen  so  that  the  right  side  of  the  first  equation 
in  (8)  vanishes.  Similarly  when  &  =  a  +  1,  then  Sa,a  must  be  chosen 
so  that  the  right  side  of  the  second  equation  in  (8)  vanishes.  By 
hypothesis  integer  (j  =  1,  •  •  • ,  m).  Therefore 

a  —  i 0  +  ix~lt*j  =^0  (j  =  l,  •••,») 

and  all  the  coefficients  of  R  can  be  computed  recursively  from  (8). 
All  the  coefficients  of  P  and  Q  can  be  computed  recursively,  except 
Pfi+i.t  and  Qa,a+l  (a,  j0  =  1, 2,  •••)  which  can  be  chosen  arbitrarily. 
Let 

PiP,  9)  =  i.fpt+W 
QiP,  q)  =  £Q«,«+iPV+l 

o=l 

by  any  given  convergent  power  series  (for  example  P,  Q  =  0),  and 
define 

—  P #  +  «.# 

Qa,*+ 1  =  Q...+*  (o,  &  ~  1, 2,  •  •  •) . 

Now,  from  (8)  and  (9)  the  coefficients  of  P,  Q,  R,  a,  S  can  be  computed 
recursively. 

Next  we  want  to  show  that  a  and  3  are  equal  as  formal  power 
series.  Since  H  is  an  integral  for  (1), 

3(p,  q)  =  H{p  +  P,  q  +  Q,  R) 

is  a  (formal)  integral  for  system  (3)  restricted  to  the  manifold  given 
by  r  =  0;  namely  for  the  system 

i>  =  i\(l  +  a)p 

Q  =  —ix(l+  S)q  . 
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Therefore,  operating  only  in  the  sense  of  formal  power  series,  we 
obtain 

4-3  =  8piX(l  +  a)p  -  8,ix(  1  +  S)q  =  0  , 
at 

or  equivalently 

(10)  p8p  -  qB ,  =  ~ap8p  +  Sq8, . 

Let 

8(p,  q)  =  2  3aJp*q?  . 

*+nt 

Equating  coefficients  in  (10),  we  have 

(11)  (a  -  P)8a,t  =  [-ap8,  +  &q8,]a.> . 

Since  a  =  a(ptf)  and  3  =  cL{pq)  are  power  series  in  the  product  pq,  it 
follows  from  (11)  by  induction  that  8tttP  =  0  for  all  a  ^  13.  There¬ 
fore  3  is  a  (formal)  power  series  in  the  product  pq,  and  this  implies 

p8p  =  qB ,  . 

Hence  from  (10) 

(a  —  S)p8,  =  0  . 

Since 

pBp  =  iXpq  +  •  •  • 

is  a  nontrivial  power  series,  we  conclude  that  as  formal  power  series 

a  {pq)  =  S(pq)  . 

Equation  (6)  can  now  be  written 

pP,  -qP,~  P  =  -a p-  a [pP,  -  qP ,}  -  iX~lX(p  +  P,  q  +  Q,  R) 
pQp  -  qQ'  +  Q  ~aq  -  a{pQp  -  gQ,}  -  iX~lY(p  +  P,q  +  Q,R) 

(  '  pR,  -  qR,  +  iX~lB?R 

=  -iX-'B'R  -  a[pRp  —  qR%)  —i\~'Z(p  +  P,q  +  Q,  R) . 

Let  X\  Y’ ,  Z’,  P’,  Q\  a '  be  the  smallest  majorants  of 

i\-lX,  ix~l  Y,  iX~lZ ,  P,  Q,  a 

respectively.  Thus 

X',  Y\  Z'  e  r(x,  y ,  z) ;  P',  Qf  €  r(p,  q)  . 

The  power  series  a'  must  still  be  shown  to  be  convergent.  Define 


TWO-DIMENSIONAL  SUBCENTER  MANIFOLDS 


341 


(13) 


P'(P,  9)  =  E(l «  ~  P I  +  8.,)  |  Pa,t  |  pV 
Q'ip,  g)  =  E(l  a  -  /3 1  +  Saf)  I  Qatf  I 

R'iip,  g)  =  E(l  a  -  fi  I  +  *«)  I  I  pV 

R  =  (R, 


(j  =  i,  •••,*») 


where  is  the  Kronecker  delta.  Clearly 

P  <  P';  0  <  O';  P«P'. 

Choose  K  >  0  such  that  for  all  nonnegative  integers  a  and  $  (except 
as  noted) 


(14) 


\a~/3-  l|-‘(|a-£|  +  8m,)£K  (a  - /3  -  1*0) 
\a-/3  +  l  j_1(| a  —  /3\  +  8af)  £  K  (a  -  0  +  1  *  0) 
\a  ~  0  +  |-1(|  a  —  0\  +  8af)  £K  (j  =  1,  •  m) . 


Prom  (12),  (13),  (14)  it  follows  that 


(15) 


P'  <  K[a'P'  +  X’(p  +  P',q  +  O',  R)  1  +  P' 

O'  «  JT[a'Q'  +  Y'(p  +  P',q  +  O',  P')]  +  G’ 

R  «  JC[|  \~l  |  BlR  +  a'R  +  Z\p  +  P'fq  +  O',  R)] 
a'p  <  a'P'  +  X'(p  +  P’,q  +  O',  S') . 


Since  P',  O',  P',  o'  all  have  real,  nonnegative  coefficients,  it  is  sufficient 
to  prove  convergence  when  p  =  q  =  £.  Let 

P'(C,  0  =  CP"(C);  0'(C,  0  =  CO"(C);  P'(C,  0  =  CP"(C); 

P'(C,  0  =  CP"(C);  Q'(C,  0  =  CO"(C);  o'(C*)  =  o"(C) . 

Prom  (15)  we  have 

P"  <  P[a"P"  +  t-lX'(C  +  CP",  C  +  CO",  CP")]  +  J?" 

0"  <  p[a''0"  +  c-‘y'(c  +  cp",  c  +  co",  cp")]  +  o" 

P"  <  K\\ X-‘  I  PlP"  +  a"P"  +  C-‘Z'(C  +  CP",  C  +  CO",  CP")] 
a"  <  a"P"  +  C"l*'(C  +  CP",  C  +  CO",  CP") . 


(16) 


Define  formal  power  series  P*  =  P*(C),  P*(0)  =  0,  etc.,  by  the  func¬ 
tional  equations 


(17) 


P*  =  K[a'P*  +  C^-HC  +  CP*,  C  +  CO*,  CP*)]  +  P" 

0*  =  P[o*0*  +  C-lP'(C  +  CP*,  C  +  CO*,  CP*)]  +  0" 

P*  =  P(|  |  P‘P*  +  o*22*  +  C"'Z'(C  +  CP*,  C  +  CO*,  CP*)] 


or 


=  o*P*  +  C-lA-'(C  +  CP*,  C  +  CO*,  CP*) . 


Comparing  (17)  with  (16)  it  follows  that 
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P"  <  P*;  Q"  <Q*;R”<R*;  a "  <  a*  . 

But  the  convergence  of  P*,Q*,R*,a*  in  (17)  follows  from  the 
implicit  function  theorem  provided  <5  in  B1  is  chosen  sufficiently 
small  to  insure  that 


Thus  from 

we  conclude  that 

To  show  that 


<  1 . 

P*,Q*,R*,a*er(  C) 


P,  Q,  R,  a  e  P*(p,  q) . 
P,Q,fter*(p,  q,r) 


write  equation  (5)  as 

P  —  — iXap  +  ixP  +  X(p  ■+•  Pf  q  +  Q,  r  +  R) 
-P,{iM  1  +  a)p  +  P)~  P,{-iX(  1  +  a)q  +  Q) 

Q  =  iXaq  -  ixQ  +  Yip  -f  P,q  +  Q,r  +  R) 

'  °  -Q,{iX(l  +  a)p  +  P)  —  Q'{-iX(l  +  a)q  +  Q} 

ft  =  BR  +  Zip  +  P,  q  +  Q,  r  +  R) 

-R„{iX(l  +  a)p  +  P}  —  R'{-iX(  1  -f  a)q  +  0} . 


Since  P,  Q,  R,  a  are  known  analytic  functions  of  (p,  q),  we  can 
regard  (18)  as  an  equation  in  the  variables  p,  q,  r,  P,  Q,  ft  and  use 
the  implicit  function  Theorem  to  solve  for  P  =  P(p,  q,  r),  etc.  Clearly 


P,Q,ft  =  0  when  r  =  0  . 


This  completes  the  proof  of  Theorem  1. 


4.  Real  analytic  invariant  manifolds.  For  real  analytic  systems 
of  equations  with  an  integral,  we  want  to  show  the  existence  of 
real  analytic  invariant  manifolds.  Theorem  1  will  be  the  prime  tool 
for  this  program. 

Consider  the  real  analytic  system  of  ordinary  differential  equa¬ 
tions  with  an  integral  H 

H{x,  y,  z)  =  -i-Mx*  +  if)  +  F{z)  +  G{x,  y,  z) 

is  -Xy  +  X(x,  y,  z) 

V  =  -Xx  +  Y(x,  y,  z) 
z  =  Bz  +  Z(x,  y,  z) 


(19) 
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where  x  and  y  are  real  scalar  variables;  z  is  a  real  m- vector  variable; 
X  is  a  nonzero  real  number;  B  is  a  real,  constant  matrix  with  eigen¬ 
values  ftit  fim ;  X,  Y,Z&  r*(x,  y,  z)  are  real  analytic ;  dim  X  = 
dim x,  etc.;  F  is  a  real  quadratic  form  in  the  components  of  z; 
G  €  r*(x,  y,  z)  is  real  analytic ;  if  is  a  real  analytic  integral  for  (1). 

Theorem  2.  If  ix-'ftj  integer  (j  =  1,  •  •  • ,  m),  then  for  system 

(19)  there  exists  a  unique ,  local,  real  analytic,  two-dimensional, 
invariant  manifold 

Ml  =  {(*,  y,  z)  1 1  x  |  +  |  y  \  <  8,  z  =  w{x,  y)} 

where  w  e  r*(x,  y)  is  real  analytic  and  is  directly  computable  ( thus 
it  is  not  necessary  to  reduce  the  differential  equations  to  a  normal 
form  such  as  (38)  below  in  order  to  compute  M1) ;  moreover,  the 
manifold  M1  is  composed  of  a  nested,  one-parameter  (p  ^  0)  family 
of  periodic  orbits,  and  as  p—* 0  the  corresponding  periodic  orbit 
goes  to  the  origin  and  its  period  goes  to  2ot  |  X-1 1. 

Proof.  The  second  part  of  this  theorem  is  relatively  easy  to 
prove.  Suppose  the  real  analytic  function  w  which  defines  Ml  is 
known.  Then 

S(x,  V )  =  -I’Ms*  +  V*)  +  F(w(x,'y))  +  G(x,  y,  w(x,  y)) 
is  a  real  analytic  integral  for  (1)  restricted  to  Ml.  Introduce 

X  =  p(l  +  b)  COB  0 

y  =  jO(l  +  6)  sin  0 

into  the  equation 

#(*,»)  =  Axp* 

to  obtain 

(20)  6  +  Aft*  +  p-*X-\P  +  6)  =  0 

A 

where 

P  =  P{q,  p,  b)  =  F(w(p(l  +  b)  cos  Q,  p(l  +  6)  sin  0)) 

&  =  Q(0,  p,  b)  —  G(p(  1  +  b)  cos  $,•••). 

We  observe  that  p-'X-'iP  +  &)  is  well  defined  for  p  =  0.  From  (20) 
by  means  of  the  implicit  function  theorem  we  can  solve  for  b  =  b(0,  p) 
with  6  having  period  2^  in  0  and  b(6, 0)  =  0,  so  that 
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*  =  (0(1  +  6(0 ,  (0 ))  COB  0 
y  =  p(l  +  6(0,  (0))  sin  0 

z  =  w(p(l  +  b(8,  p ))  cos  0,  (0(1  +  6(0,  (0))  sin  0) 

is  a  representation  of  the  one-parameter  (p  0)  family  of  nested 
periodic  orbits  which  comprise  M J. 

Continuing  in  this  manner,  it  is  not  difficult  to  prove  the  asser¬ 
tion  concerning  the  periods  of  these  periodic  orbits.  See  [3]  for  the 
details. 

For  a  more  complete  analysis  of  the  differential  equations  (19) 
restricted  to  M1,  one  can  show  the  existence  of  a  real  analytic  change 
of  variables 

x-p  +  P(p,  q) 

V  =  q  +  Q(P ,  q) 
z  =  r  4-  R(p,  q) 

such  that  on  Ml  (given  by  z  =  R(p,  q)) 

P  =  Ml  +  a(P*  +  ?*))? 

$  =  -Ml  +  a(P*  +  fl*))P 

where  P,  Q,  R,  a  e  P*( p,  q)  are  real  analytic ;  a  is  a  real  power  series 
in  the  quadratic  form  p*  +  q*.  In  this  formulation  the  periodicity  of 
the  solutions  of  (21)  is  obvious  (the  circles  p*  +  q*  =  (0*  in  the 
(p,  q)-plane  are  orbits)  and  the  assertion  concerning  the  periods  of 
these  solutions  is  immediate.  The  details  of  this  analysis  are  given 
in  §5  below. 

We  now  prove  the  first  part  of  Theorem  2.  Primes  for  the  rest 
of  this  section  will  denote  variables  and  functions  in  a  complex 
setting  as  opposed  to  unprimed  variables  and  functions  in  a  real 
setting.  Introduce  the  change  of  variables 

x  =  +  v') 

(22)  y  =  ^~(-*'-»V) 

z  —  Bz' 

so  that  system  (19)  transforms  to 

H\x’,  y\  z')  =  ixx'y '  +  F\z')  +  G(x’,  y\  z') 
i’  =  ixx'  +  X'(x',  y',  z’) 
y'  =  -i\y’  +  Y'(xf,  y\  zf) 
z'  =  B'z!  +  Z'( x\  i/*, «') 


(23) 
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where 

H'tf,  z')  =  +  j/0,  -  »V),  Bz')  ; 

J5'  =  B-'BB  is  in  Jordan  canonical  form;  X’,  Y',  Z’  are  related  to 
X,  Y,  Z  in  the  usual  way.  Theorem  1  asserts  the  existence  of  an 
invariant  manifold 


(24)  ‘  =  {(X'’ V>' Z>)  1  X’  =  P'  +  P'(P’’  ?'>, 

*  y'  =  q'  +  Q’(p',  90,  z'  =  R'{P\  90} 

where  P',  Q',  R'  e  r*(p',  q').  But  the  change  of  variables 

x'  =  p'  +  P'(p',  q') 

V'  =  q'  +  Q'(p\  q') 

can  be  inverted  to  give 

p’  =  p’{x',  y')  =  x'  +  •  •  • 
q'  =  q'(xr ,  y')  =  y'+  •  •  • 

so  that 

(26)  M'*  =  {(x',  y,,zf)\\xf\  +  \y'l  <d,  z>  =  y')) 

where 

w'(*',  VO  =  R'{p'{*,  V'),  q'{x’,  I/O)  , 

which  implies  w'  e  r*(x',  y').  With  the  representation  given  by  (24) 
it  is  not  clear  that  Mn  is  unique,  but  with  the  representation  (25) 
we  can  easily  show  uniqueness  for  Mn.  Since  the  manifold  is 
invariant, 

W  =  B'w’  +  Z'(x',  y',  w  0  , 

equivalently 


(26) 


iXx'w —  iXy’w’,.  —  B'w' 

=  Z'(x',  i/O  w')  -  w'„X'( x',  y’,  w’)  -  w'rY'(x',  y't  w’)  . 


Expanding  to'  as  a  power  series  and  equating  coefficients  in  (26),  we 
have 


2?)  (iX(a  -  0)1  -  B'W',' 

=  [Z'(x',  irO  -  *..X'(x',  i/,  w')  -  w'r  Y'itf,  y', 

where  I  is  the  identity  matrix.  Since  by  hypothesis  integer 

(j  =  1,  •*., m)  where  (tu  •••,[*„  are  the  eigenvalues  of  B*,  it  follows 
that 
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det  (i\(a  -  0)1  —  B')  ^  0  . 

Since  the  right  side  of  (27)  depends  only  on  the  coefficients  of  w' 
having  order  less  than  a  +  0,  we  conclude  that  w'  is  uniquely 
determined  by  (27)  and  that  the  manifold  Mn  is  unique.  Define 

(28)  w(x,  y)  =  Bw'(^—{-ix  -  y),  ^L(x  +  iy))  . 

Then 

(29)  tb  =  Bw  +  Z(x,  y,  w)  , 
equivalently 

(30)  XVWx  ~  XXW>  ~  BW 

!  =  Z(x,  y,  w)  —  wmX(x,  y,  w )  -  wtY(x,  y,  w)  . 

The  function  w  given  by  (28)  satisfies  (30).  Suppose  w*er*(x,  y)  is 
another  solution  to  (30).  Then  w*  determines  a  function  w*'  e  /*( x',  y') 
given  by 

(31)  w*'(x',  y')  =  B-lw*(J^(ix'  +  y'),  JL(-x'  -  iy'))  . 

But  since  w*  satisfies  (30),  it  also  satisfies  (29).  Then  from  (31)  we 
conclude  that  w*'  satisfies  (26),  hence  w *'  =  w'  and  because  (28)  and 

(31)  are  inverse  relations  to  each  other,  w*  =  w.  Therefore  we  con¬ 
clude  that  (30)  has  a  unique  solution  belonging  to  the  class  r*(x,  y). 

It  remains  to  be  shown  that  mi  is  a  real  analytic  function.  If 
we  expand  the  left  side  of  (30)  in  a  power  series,  we  obtain 

Xyw,  —  Xxw,  —  Bw 

(32)  =  £\aw<tl/,xB-y+1  -  ^lX0wa,fxa+ly?-1  -  Y,Bwa,txty? 

=  E(M«  +  l)w«+i.#-»  -  MB  +  1  )wa-t,f+l  -  Bw.Jx'yt . 

For  any  power  series 

/ = 

define 

/.  =  [/!•  =  IE/«,>*V]«  =  (/..0,  /*-».»,  *  ’  *  ,  /<>,•) 

for  all  nonnegative  integers  n.  Thus  if  /  is  a  vector,  then  dim/.  = 
(n  +  l)dim/.  If  /  and  f  are  two  power  series  related  by 

/'<*',  1/)  =  fix,  y)  =  /(^< ix’  +  y'),  ~^(-x'  -  iy ')) 

f(x,  y)  =  fix’,  y1)  =  /'( Jy(-ix  -  y),  ^J-(x  +  iy))  , 
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then  for  each  nonnegative  integer  n  there  exists  a  matrix  L%, 
det  L„  0,  such  that 


=  LJn  . 

Clearly  the  matrix  L%  depends  only  on  dim  /;  other  than  this  L%  is 
independent  of  f. 

Now,  from  (32)  we  obtain 


(33)  [Xywx  —  Xxw,  —  Bw]%  =  B,wk  (n  =  2,  3,  •  •  •) 

where 

Bh  =  super-diag (— X7,  —  2 X7,  •  •  •,  —nxl) 

+  diag  (—7?,  -B,  •  •*,  -B) 

+  sub-diag  (n\7,  (n  —  1)X7,  •  •  • ,  x7) . 

with  7  =  identity  matrix,  dim  7  —  dim  B.  In  a  similar  fashion  we 
obtain 


[iXx'w'a.  —  iXy’w'y,  —  B'w']n  =  Bfnw'n  (n  —  2,  3,  •  •  •) 


where 

B:  =  diag  (inXl~B\  i(n  -  2)XI-B\  i(n  -  4)Xl-B',  ••, -inXl-B') . 
Since  det  (ikxl  —  B')  =£  0  for  any  integer  k,  it  follows  that 

det  B'%^Q  (n  =  2, 3,  •  •  •)  • 

A  straight  forward  calculation  shows  that 


.  3  N  3  ..  ,  3  ...  3 

Xy— - Xx—  —  iXx  -r-—  —  iXy  -  ■  , 

ox  oy  ox  oy* 

and  from  this  fact  along  with  (28)  we  conclude 


Xyw,  —  Xxw,  —  Bw  =  B{iXx'w'a,  —  iXy'w',,  —  B’w] 

(34)  [Xyw,  —  Xxw,  —  Bw]n  =  La1[B{\Xx'wa,  —  iXy'w',,  —  B'w')]n 
B,wn  =  L-'B.B’.w'.  =  L~lB%ffnLnw%  , 


where 

Bn  =  diag  (B,  B,  •  •  • ,  B) . 


Since  (34)  can  be  shown  valid  for  all  wn  e  it  follows  that 

B.  =  Lt'B.B'.L. 
det  Bm  =  det  B » •  det  B".  ^  0  . 


(85) 
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From  (30)  and  (38)  we  obtain 

(36)  B%w%  =  [Z(x,  y,  w)  -  wMX(x,  y,  w)  -  w,Y(x,  y,  to)]. 

(n  =  2, 3,  •  •  •) . 

Finally,  from  (35)  and  (36)  we  conclude  the  coefficients  of  w  are 
determined  recursively  from  (30),  or  equivalently  (36),  and  since 
B%  (n  =  2, 3,  •  •  •)  is  a  real  matrix,  B^1  is  also  real  and  w  e  /'*(*,  y)  is 
real  analytic.  Thus  we  have  shown  that  the  function  w  which 
defines  the  manifold 

M 1  =  {(*,  y,  z)  |  j  x  |  +  |  y  |  <  8,  z-  w{x,  y)} 

is  uniquely  and  recursively  determined  from  (30),  (35),  (36)  as  a  real 
power  series.  The  fact  that  w  e  /*( x,  y)  is  guaranteed  by  the 
existence  of  the  real  analytic  integral  H.  This  completes  the  proof 
of  Theorem  2. 

5.  Real  analytic  changes  of  variables.  In  this  section  for  a 
real  analytic  system  of  equations  we  construct  a  real  analytic  change 
of  variables. 

Theorem  3.  If  &  integer  (j  =  1,  •  •  •,  m)  then  for  system 

(19)  there  exists  a  real  analytic  change  of  variables  ( not  unique) 


x  =  p  +  P(p,  q ) 

v  =  q  +  Q(p,  q) 

z  =  r  +  R(p,  q) 

i>  =  \(1  +  a(p*  +  q'))q  +  P(p,  q,  r) 

-Ml  +  a(P *  +  q*))p  +  Q(p,  q,  r) 
r  =  Br  +  R{p,  q,  r) 

P,  Q,  R  =  0  when  r  —  0  , 

where  a  =  a(p*  +  ql)  is  a  real  analytic  function  of  the  quadratic 
form  px  +  q'  with  no  constant  term  in  its  expansion  about  the 
origin;  P,  Q,  Re r*(p,  q)  are  real  analytic;  P,£},  Be  r*(p,  q,  r)  are 
real  analytic. 

Proof.  By  means  of  the  change  of  variables  (22)  we  obtain 
system  (23).  From  Theorem  1  we  can  find  a  change  of  variables 


(37) 

such  that 

(34) 
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*  =  p'  +  P'(p',  q') 

(39) 

v'  =  q'  +  Q'(p',q') 

*  =  r'  +  Rtf,  q ') 

such  that 

V  =  *M1  +  a’{p’q'))p'  +  P(p',  q>,  r') 

(40) 

4'  =  -i\(  1  4-  a'(p'9'))9'  +  0'(p',  9',  tO 
P  =  PV  +  #(p',  g\  P) 

P',  O',  &'  =  0  when  r'  =  0 


where  a'  =  a'(p'g')  is  a  convergent  power  series  in  the  product  p'q' 
with  no  constant  term;  P',  O',  1?  e  P*(p',  q') ;  P,  0,  8 e  P*(p',  g',  r'). 
In  the  construction  of  the  change  of  variables  (39)  one  must  restrict 
the  choice  of  power  series  P',  Q'  (see  (8)  and  (9))  to  be  transforma¬ 
tions  via  (22)  of  real  power  series.  We  will  make  this  condition 
explicit  in  (43)  below.  Define 

P  =  '^2^'  "*■  90 

(41)  q  =  -  iq') 

r  —  Sr1 . 

Then  from  (22),  (39),  (41)  we  obtain  the  change  of  variables  (37) 
where 


0^  +  9*)  =  i(p*  +  g*)) 

P(p,  q)  =  ^%P'{-^-ip  -  q),  _^(p  +  iq)j 

(Up,  9)  =  -~^P'{^i-ip  -  q),  Jj(p  +  iq)) 

~vti<r(v^~ip  -  q)’  VT(P  +  *>) 

R(p,  9)  =  SR'(^j(-ip  -  q),  ^<P  +  »9))  . 


In  particular  we  require  P,  Q  given  by 
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P(p,  Q)  -  q),  ^-(p  +  ig)) 

+ V2®(v^~iP  ~  q)’  V2iP  +  ^}) 

(43)  -  l  */  1  1  \ 

QiP,  Q)  -  -^=P'(^-(-tp  -  q),  y-{p  +  iq)) 

-V^7*~**  "  9)*  V2(P  +  *>) 

to  be  real  analytic  functions  of  (p,  q).  The  functions  P,  Q,  R  are 
given  by  formulas  similar  to  those  for  P,  Q,  R  in  (42). 

It  remains  to  show  that  a,  P,  Q,  R,  P,  Q,  R  are  real  analytic 
functions.  From  (19),  (37),  (38)  we  obtain 

qP,  -  pP,  -  Q  =  -aq  -  a{qP,  -  pP,)  +  X~lX(p  +  P,q  +  Q,R) 

(44)  qQ,  -  pQt  +  P  =  ap  -  a{qQ,  -  pQ,}  +  X-‘F(p  +  P,  q  +  Q,  R) 
qRr  —  pRt  —  X~lBR  -  - a{qR ,  -  pR,}  +  \~lZ(p  +  P,  q  +  Q,  R) . 

Following  a  procedure  analogous  to  the  computation  of  w  in  (30)  we 
can  compute  the  coefficients  of  a,  P,  Q,  R  recursively  from  (44).  Since 
X,  Y,  Z  are  real  analytic  and  X  is  a  real,  nonzero  scalar,  it  follows 
that  a,  P,  Q,  R  are  real  analytic.  This  completes  the  proof  of 
Theorem  3. 
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